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The Fenna-Matthews-Olson (FMO) antennae complex, responsible for light harvesting in green 
sulfur bacteria, consists of three monomers, each with seven chromophores. Here we show that 
multiple subsystems of the seven chromophores can transfer energy from either chromophore 1 or 6 
to the reaction center with an efficiency matching or in many cases exceeding that of the full seven 
chromophore system. In the FMO complex these functional subsystems support multiple quantum 
pathways for efficient energy transfer that provide a built-in quantum redundancy. There are many 
instances of redundancy in nature, providing reliability and protection, and in photosynthetic light 
harvesting this quantum redundancy provides protection against the temporary or permanent loss 
of one or more chromophores. The complete characterization of functional subsystems within the 
FMO complex offers a detailed map of the energy flow within the FMO complex, which has potential 
applications to the design of more efficient photovoltaic devices. 
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LETTER 

Most photosynthetic organisms have light-harvesting 
antennae complexes composed of pigment molecules, 
called chromophores, whose electrons are excited by sun- 
light [ij. The excitation energy is transferred from the 
antennae complexes to reaction centers with nearly a 
100% quantum efficiency Q. For many years this trans- 
fer was thought to occur by a classical hopping mecha- 
nism, in which the excitation energy transfers from one 
chromophore to another in a downhill fashion until it 
reaches the reaction center 0, 3 , but recently new ex- 
perimental evidence suggests that this transfer occurs 



by a quantum coherent mechanism [f| . Quantum co- 
herence in the Fenna-Matthews-Olson (FMO) antennae 
complex of green sulfur bacteria has been shown by two- 
dimensional Fourier-transform electronic spectroscopy to 
last over 500 fs at ambient temperatures and over 300 fs 
at physiological temperatures [7H9|. 

In this Letter we show that photosynthetic light 
harvesting exhibits quantum redundancy. By quan- 
tum redundancy we mean that photosynthetic anten- 
nae complexes have multiple efficient quantum path- 
ways for transferring energy to the reaction center. 
These pathways involve different subsets of the an- 
tennae's chromophores, which form functional subsys- 
tems. The antennae complex known as the FMO com- 
plex is a trimer where each monomer consists of seven 
bacteriochlorophyll-a chromophores [10] . Previous re- 
search has shown that in each monomer initial excitation 
can occur on either chromophore 1 or chromophore 6 and 
that this energy must move to chromophore 3, which is 
most closely coupled with the reaction center Here 
we show that within each monomer of the FMO complex 
there exist many subsets of the seven chromophores with 
efficiencies close to or even better than the efficiency of 
the entire set of seven chromophores. Quantum redun- 
dancy contributes a robustness to nature's quantum de- 
vice with potential survival benefits. 

The computations presented here reveal that the func- 
tional subsystems achieve their efficiencies by a quantum 
mechanism similar to that of the full system including 
the roles of entanglement and environmental noise. We 
assess each subsystem's entanglement by a global entan- 
glement measure 12 , 13| based on the squared Euclidean 
distance 14-1^. The study of these chromophore sub- 
systems gives more information about the role of each 
chromophore in the energy transfer in the whole FMO 
complex 23-3l|. This more in-depth understanding can 
provide insight to other antennae complexes and ulti- 
mately be applied to create synthetic solar cells that rival 
the efficiency of nature. 
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We consider a single monomer of the FMO complex 
with M chromophore sites with the Hamiltonian 

M 

H = J2 huj 3 <7 t a 7 + £ nv o,l{ u J a l + a j a f)> (!) 

J'=l J# 

where crT {vj) creates (annihilates) a single excitation 
on chromophore j. The site energy of each chromophore 
is hujj while the coupling constant between the pair of 
chromophores j and I is Hvjj. 

The time evolution of the system's density matrix is 
governed by the quantum Liouville equation 



±D = - i -[H,D] + L { D) 
where D is the density matrix 

it./ 



(2) 



(3) 



pf are the elements of the density matrix D in the basis 
set of wavefunctions {^fc}, and L is the Lindblad op- 
erator, which accounts for interactions of the M chro- 
mophores with the environment. At t = we initialize 
the density matrix with a single excitation (exciton) on 
either site 1 or 6. Because the Hamiltonian does not mix 
singly excited wavefunctions with other wavefunctions, 
the set {^k} only includes M wavefunctions formed from 
sequentially considering a single excitation on each of the 
M chromophores. 

The Lindblad operator 3(1 32| is the sum of three op- 
erators that account for dephasing, dissipation, and loss 
to the reaction center (sink) 



L(D) = L deph (D) + L diss (D) + L sink {D) (4) 



where 



L deph (D)=aJ2 2 (k\D\k)\k)(k\-{\k)(k\,D}, (5) 

k 

L diss (D) = W\D\k)\g){9\ - {\k){k\,D}, (6) 

k 

L sin k(£')=27(3|£>|3)| S )( S |-7{|3)<3|,£>}. (7) 

The \g) denotes the state in which each of the M chro- 
mophores is in its ground state, the \k) represent the 
M excited states with each state having a single excita- 
tion on one of the M chromophores, the state \s) denotes 
the reaction center (sink), and |3) indicates the excited 
state of the third chromophore multiplied by the ground 
states of the other M — 1 chromophores. The particu- 
lar form of the Lindblad operators chosen here follows 
the operators selected in previous studies of the FMO 
complex [13]. Further refinements of the system's inter- 
action with the environment will not substantively af- 
fect the general conclusions made here about the role of 



functional subsystems within the FMO complex. The 
model incorporates two types of noise that affect the en- 
ergy transfer. First, dissipation transfers energy from 
the chromophores back into the environment, wasting 
that energy. Second, dephasing dampens the coherence 
within the chromophore system. Finally, the Lindblad 
operator L s ; n k(-D) is responsible for transferring energy 
from the third chromophore into reaction center, which 
acts as an energy sink and becomes a M + 1 site in the 
density matrix. The rate parameters a, /?, and 7 were 
chosen in atomic units to be 1.52xl0 -4 , 7.26xl0 -5 , and 
1.21xl0~ 8 , respectively. These values are consistent with 
those employed in Ref. [30]. 

Entanglement is a correlation that cannot occur in a 
classical system; Bell [HI defined entanglement as "a 
correlation that is stronger than any classical correla- 
tion." Parts of a quantum system become entangled 
when the total density matrix for the system cannot be 
expressed as a product of the density matrices for the 



parts [H H 13. Global entanglement [12J, [lj, |34| 

the system is measured by the squared Euclidean distance 
between the density matrix p and the nearest classical 
density matrix £ [l~4T - [l9j 



(8) 



In the basis of M wavefunctions, in which each wave- 
function has a single excitation on one of the M chro- 
mophores, this measure corresponds to the sum of the 
squares of the off-diagonal elements of p, or the sum of 
the squares of the concurrences. The squared Euclidean 
distance is zero if and only if the excitations (excitons) on 
the chromophores are not entanglement. Both squared 
Euclidean distance and quantum relative entropy [29I l3~il] . 
another measure of global entanglement, are Bergmann 
distances. 

In Figs. 1-3 and Table I we study subsystems of 
the FMO complex, in which certain of the seven chro- 
mophores are "turned off" and the energy transfer is 
tracked within a smaller set of chromophores. For the 
full system of seven chromophores we employ the 7x7 
Hamiltonian given in Ref. 30]; chromophores are re- 
moved from the full system by deleting rows and columns 
of the Hamiltonian matrix. We denote the subsystem by 
the numbers of the chromophores retained; for example, 
the subsystem with chromophores 1,2, and 3 is denoted 
as 123. Figure 1 shows the exciton populations of chro- 
mophores 1, 2, 3, and the reaction center for reduced 
system 123 and the full system. We observe that the 
two systems are quite similar, and yet they exhibit some 
important differences. The reduced system features a 
slightly slower decay of the exciton populations in chro- 
mophores 1 and 2, and most significantly, it transfers 
energy from chromophore 1 to the reaction center more 
quickly than the full system. Because chromophores 4-7 
draw some of the excitation energy in the full system, the 
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FIG. 1. Populations of chromophores 1-3 and the reaction 
center for (A) the subsystem 123 and (B) the full FMO sys- 
tem 1234567. Subsystem 123 of the FMO complex transfers 
energy more efficiently from chromophore 1 to the reaction 
center (sink) than the full FMO system. 



populations of chromophores 1-3 and the reaction center 
are lower in the full system than in the reduced system. 

The efficiency of energy transfer from site 1 to the 
sink is compared in Fig. 2 for several subsystems, 123, 
1234, 12345 (not shown because its efficiency is similar 
to 1234), and 123456, as well as the full system. We 
observe that subsystems 123, 1234, and 12345 are more 
efficient than the full system even though efficiency de- 
creases as we add chromophores from 123 to 123456. The 
addition of chromophores 4, 5, and 6 increases the num- 
ber of sites for the energy to enter and thereby decreases 
the population of every site, including the reaction cen- 
ter. The addition of chromophore 7 to 123456, however, 
improves efficiency; because chromophore 7 has a lower 
energy than chromophores 5 and 6, it offers a new down- 
hill quantum path through chromophore 4, even lower in 
energy, to chromophore 3 and the reaction center. The 
importance of chromophore 7 can also be seen from Ta- 
ble I, which shows that the subsystem containing chro- 
mophores 12347 is the most efficient system, again be- 
cause of the downhill path created from chromophore 7 
to the reaction center. 
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FIG. 2. Population in the reaction center for several chro- 
mophore subsystems. The subsystems 123, 1234, and 12345 
(not shown) are more efficient than the full FMO system at 
transferring energy from chromophore 1 to the reaction center 
(sink). 



Table I reports the efficiencies of subsystems with 
the initial excitation on either chromophore 1 or chro- 
mophore 6. Subsystems with both chromophores 1 and 
6 are included in both categories. The table shows that 
most of these subsystems are more efficient than the full 
system at transferring energy from either chromophore 1 
or chromophore 6 to the reaction center. However, only 
one subsystem 123467 is more efficient than the full sys- 
tem when efficiency is assessed with equal weight given 
to the excitation pathways from chromophores 1 and 6 
to the reaction center. Although this model suggests 
that chromophore 5 is superfluous, it may have other 
necessary duties in nature that are not considered in 
the present model. The superfluity of chromophore 5 
is predictable from the Hamiltonian of the system, as 
the site energy for chromophore 5 is much higher than 



the other chromophores 11] . However, it is worth not 



ing that this is only the case in the FMO complex for 
Prosthecochloris aestuarii , which we are studying in 
this Letter, while in the FMO complex for Chlorobium 
tepidum, chromophore 5 is not the highest in energy and 
its exclusion would most likely make the system less ef- 
ficient. Very recent crystallographic and quantum chem- 
istry studies 35, 36} indicate that there is likely an eighth 
chromophore in the FMO of Prosthecochloris aestuarii, 
which due to sample preparation is not present in the 
ultrafast spectroscopic studies An eighth chro- 

mophore would further support the existence and im- 
portance of functional subsystems with the seven chro- 
mophores in recent ultrafast experiments serving as a 
functional subsystem of the full system in nature. 

The subsystems, we observe, function by a quantum 
mechanism similar to that observed in the full system. 
Two features of this mechanism, which have received sig- 
nificant attention for the full system in the literature |20l — 
3l| . are: (i) dephasing from environmental noise and (ii) 
entanglement of excitons. As in the full system, the 
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TABLE I. Populations of the reaction center of different chro- 
mophore subsystems at 2 ps. 

Reaction Center Population 
Initial Excitation System No Dephasing Dcphasing 
Chromophore 1 123 0387 0592 



1234 0.161 0.509 

12345 0.144 0.506 

123456 0.120 0.445 

1234567 0.459 0.498 

12347 0.347 0.564 

123467 0.314 0.514 

Chromophore 6 346 0.010 0.090 

367 0.009 0.026 

3467 0.048 0.486 

3456 0.045 0.305 

34567 0.103 0.445 

234567 0.086 0.434 

1234567 0.139 0.433 

123467 0.074 0.456 




Time (fs) 



FIG. 3. Entanglement of excitons for several chromophore 
subsystems as a function of time. A measure of global entan- 
glement is shown as a function of time (fs) for several chro- 
mophore subsystems, 123, 1234, and 12346, as well as the full 
system 1234567. Like the full system, the subsystems exhibit 
entanglement of excitons between the chromophores. 



dephasing assists transport of the excitons from chro- 
mophore 1 or 6 to the reaction center by facilitating 
the distribution of the exciton population from its ini- 
tial chromophore. The efficiencies of different subsystems 
with and without dephasing are reported in Table I. The 
perceived enhancement is a characteristic of systems of 
heterogeneous chromophores, that is chromophores with 
different site energies. We also observe in Fig. 3, dis- 
playing the entanglement in several subsystems as well 
as the full system, that the entanglement of excitons ex- 
hibits similar behavior in the reduced systems as in the 
full system. We compute the entanglement to demon- 
strate that the transfer of energy for both reduced and 
full systems within the model occurs by the same quan- 
tum mechanism. (It is not the purpose of the present 
Letter to examine the relationship between entanglement 
and efficiency, which has been explored in previous work 
including a recent paper by one of the authors |37[-) The 
addition of each chromophore slightly increases the total 
entanglement of the system, as measured by the squared 
Euclidean distance, by adding extra off-diagonal elements 
to the density matrix. Each system shows that the en- 
tanglement lasts for a few hundred femtoseconds, which 
matches previous experimental data as well as predictions 
from other theoretical models 0, [38| . 

Efficient photosynthetic light harvesting, according to 
recent experimental evidence, likely occurs by a quan- 
tum rather than a classical mechanism. In this Let- 
ter we have identified multiple functional subsystems of 
the FM O light-harvesting complex whose efficiencies are 
comparable to the efficiency of the full FMO complex. 
Many of these subsystems are, in fact, more efficient in 
transferring energy from either chromophore 1 or 6 to the 
reaction center, and one of these subsystems is even more 
efficient in transferring energy from both chromophores 1 
and 6 to the reaction center. These functional subsystems 



provide evidence for quantum redundancy in photosyn- 
thetic light harvesting. There exists multiple quantum 
pathways or circuits that are capable of transferring en- 
ergy efficiently to the reaction center. For example, sub- 
system 123 transfers energy from chromophore 1 to the 
reaction center with 59.2% efficiency (after 2 ps), and 
subsystem 3467 transfers energy from chromophore 6 to 
the reaction center with 48.6% efficiency. This built-in 
redundancy likely provides benefits to the natural sys- 
tem because damage to any chromophore save 3 will not 
disrupt light harvesting. The characterization of func- 
tional subsystems within the FMO complex offers a de- 
tailed map of the energy flow within the FMO complex 
with potential applications to the design of more efficient 
photovoltaic devices. 
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